function check_Bm_period(CONSTS, plot_data)

    max_m = 30;

    m_vec = (-max_m:1:max_m)';
    %m_vec = (1:1:1)';
    Bm_vec = zeros(size(m_vec));
    Cm_vec = zeros(size(m_vec));
    for im = 1:size(m_vec, 1)
        m = m_vec(im);

        [coef_B, coef_C] = func_coef_B_C(m, CONSTS);

        Bm_vec(im) = coef_B; % m^2 * coef_B;
        Cm_vec(im) = coef_C;
    end

    if (plot_data)
        figure; plot(m_vec, Bm_vec, 'b.-', m_vec, Cm_vec, 'r.-')
        title('B_m(m), C_m(m)')
        xlabel('m'); ylabel('B_m'); %ylim([0 8e-13]);
        legend('B_m', 'C_m', 'location', 'SouthEast');
    end

end
